function A = getAm(N,S,alpha,psi,delta)

% ss = linspace(1,max(S),max(S))';
% aux = psi.^-delta;
% auxs= psi.^(-delta.*S);
% a = psi.^((1/(1-alpha)-delta)).^ss;
% aa = cumsum(a,1);
% ss = repmat(ss,1,size(N,2));
% sss = logical(ss==repmat(S,max(S),1));
% int = sum(aa.*sss);
% A = N.*int.*(1-aux)./(aux.*(1-auxs));

A = N.*fintA(S,alpha,psi,delta);

end
